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Abstract 

The difference in the proper azimuthal periods of revolution of two standard clocks in direct 
and retrograde orbits about a central rotating mass is proportional to J/Mc 2 , where J and 
M are, respectively, the proper angular momentum and mass of the source. In connection 
with this gravitomagnetic clock effect, we explore the possibility of using spaceborne standard 
clocks for detecting the gravitomagnetic field of the Earth. It is shown that this approach to 
the measurement of the gravitomagnetic field is, in a certain sense, theoretically equivalent to 
the Gravity Probe -B concept. 
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1 Introduction 



Currently there is considerable interest in the development of highly stable clocks for 
space applications This circumstance provides the impetus to investigate further 
a certain remarkable gravitomagnetic clock effect || within the framework of general 
relativity. 

The Newtonian theory of gravitation may be thought of as the nonrelativistic theory 
of the gravitoelectric field E 9 , while general relativity involves - among other things - 
the gravitomagnetic field B 9 as well. These notions are ultimately based on the close 
formal analogy between Newton's law of universal gravitation and Coulomb's law of 
electricity. The gravitomagnetic clock effect involves a certain characteristic temporal 
structure around rotating bodies. The elucidation of the various aspects of this effect 
is the main subject of this paper; in particular, we consider the theoretical problem of 
measuring the non-Newtonian gravitomagnetic field by means of ultra-stable clocks in 
space. 

The standard tests of the Einstein theory of gravitation can be accounted for by post- 
Newtonian gravitoelectric corrections: gravitational redshift, perihelion precession of 
Mercury, bending of light in the field of the Sun, and Shapiro's radar time delay. The 
investigation of gravitomagnetic effects in general relativity began with the work of de 
Sitter, Thirring, and Lense and Thirring ||. However, the idea of a gravitomagnetic 
field generated by mass current dates back to the last century when developments in 
electromagnetism suggested a generalization of Newton's theory of gravitation along the 
lines of electrodynamics |4|, H. Holzmiiller and Tisserand M tried to explain the 
perihelion excess of Mercury by taking into account solar gravitomagnetism. The anal- 
ogy with electrodynamics has been the subject of many studies || [| [LjJ; in fact, any 



theory that brings together Newtonian gravitation and Lorentz invariance must contain 
gravitomagnetism in some form ]TTT| . At present, general relativity is consistent with all 



observational data and within its framework the theoretical development of gravitoelec- 
tromagnetism has reached a certain level of maturity |12|. On the experimental side, 
moreover, the main mission of NASA's Gravity Probe -B will be the direct measurement 
of the gravitomagnetic field of the Earth via superconducting gyroscopes carried by a 
drag-free satellite in a polar orbit about the Earth [|TI| . 

In this paper, we choose an astronomical body such as the Earth for the sake of con- 
creteness and concentrate our attention on the theoretical possibilities for measuring its 
gravitomagnetism by means of spaceborne clocks. To this end, the exterior gravitational 
field of the Earth may be considered in the linear approximation with a metric of the 
form 



-ds 2 = -c 2 (l - 2GM ^dt 2 + ( 1 + 2GM ')5i j dx i dx j - 4 _ ( f e m J i x j dx k , (1) 
V c 2 p / V c 2 p J crpr 



where M and J are the mass and angular momentum of the body, respectively. Here 
p is the isotropic radial coordinate in an underlying Cartesian coordinate system p = 



2 



(x l ) = (x, y, z) and we choose an orientation such that J = Jz. The Lense-Thirring part 
of this metric has the form — 4c _1 (A 9 • dx)dt, where 

A9 ~ C P 3 ^ 

is the gravitomagnetic vector potential. A free ideal gyroscope held at rest at a fixed 
position p in space would precess in the gravitomagnetic field at a rate given by 



n P = — r 

c z p a 



3(J.p)p-p 2 J (3) 



in this approximation. The gravitomagnetic field is then defined byB 9 = cf2p = Vx A g , 
while the gravitoelectric field is E 9 = — V$jv, where $at = GM/p is the Newtonian 
potential, as expected. 

In the following sections 2 and 4, the Schwarzschild radial coordinate r is employed, 

r = K 1 + 2-rJ ' & 

moreover, we usually set G = 1 and c = 1 for the sake of simplicity - except where 
indicated otherwise. The gravitomagnetic clock effect under consideration in this paper 
is described in terms of azimuthal closure for simple geodesic orbits in sections 2-4. 
This gravitomagnetic effect may be thought of in terms of a certain limiting form of 
the gravitational Aharonov-Bohm effect; indeed, this connection leads to a discussion 
of holonomy in section 5. A preliminary examination of the experimental possibilities is 
provided in section 6. The measurement of the gravitomagnetic clock effect appears to 
be beyond present experimental capabilities by almost an order of magnitude. 



2 A Gravitomagnetic Clock Effect 

A standard clock by definition measures proper time along its worldline. Imagine such 
clocks in the stationary spacetime outside a charged rotating source; in fact, we take 
this region to be the exterior Kerr-Newman spacetime for the sake of simplicity. 

Consider equatorial circular geodesic orbits in this spacetime. We are interested in 
timelike orbits that are stable against radial perturbations. This is possible for circular 
orbits only beyond a certain radius. Let t + be the period of revolution in the same 
(opposite) sense as the rotation of the source for the stable circular geodesic path with 
fixed radial coordinate r in standard Schwarzschild-like coordinates. It can be shown 
that 

t± = T ± 2tt- , (5) 
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where T = 2ti/ujo, u is the modified "Keplerian" frequency given by 



u 



GM GQ' 



c 2 r 4 



(6) 



Q denotes the charge, and a = J/Mc is the Kerr parameter [14]. The coordinate time t 
is the proper time of static asymptotically inertial observers that are infinitely far from 
the source. It follows from equation (|5]) that 



47T 



J 



Mc 2 ' 



(7) 



The operational significance of this interesting result is doubtful, since light signals are 
required, for instance, to carry information about the orbit to the distant clocks. To 
ameliorate this situation, let us consider instead the proper periods, r±, of such circular 



paths measured by free orbiting standard clocks. It is possible to show that |14 

3GM 



T± 



To 



c 2 r 



+ 2 



so that 



GQ 2 

c 4 r 2 



An 



±2-u c 

c 



J 



Mc 2 1 



(8) 



(9) 



hence, an experimental determination of the left side of equation (0) would lead to the 
measurement of J/M for the source. Specifically, let us write equation @ as 



T-L 



An 



J 



Mc 2 



U(r), 



(10) 



-i 



(t + + r_)/T and U(r) approaches unity as r — > oo. Under physically 



where U 

realistic conditions, U is a monotonically decreasing function of r; in fact, for 

GM 



c 2 r 



< 1 



U (r) has an expansion of the form 

Q 2 



U(r) 



3^ /27 



8 GM 2 



135 1 c 2 J 2 



9 Q 2 



2 C.'-'.l/ 1 2 GM 2 



(12) 



The exact result is valid for the Kerr-Newman geometry. We are interested, however, 
in the exterior field of a rotating astronomical body; therefore, r GM/c 2 and hence 



An 



J 



Mc 2 



(13) 



This is a remarkable relation that could, in principle, be used to measure J/M directly 
for an astronomical body; in fact, for the Earth r + — t_ ~ 2 x 10 _7 sec, while for the Sun 
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r + — t_ ~ 10~ 5 sec. Equations (^) and (|yD indicate a general feature of time in the field 
of a rotating source; this paper is devoted to an account of this gravitomagnetic clock 
effect and its possible observational significance. 

The general result that it takes longer for a free test particle to go around a rotating 
mass in its equatorial plane in the prograde direction than in the retrograde direction 
is a remarkable fact that is in conflict with the notion that a rotating mass drags space 
around with it; in fact, such a "Machian" concept is in conflict with general relativity 
and must be abandoned |3], [y|. Moreover, the "dragging of local inertial frames" has 
often been used as a metaphor for the gravitomagnetic precession of ideal test gyro- 
scopes; however, it must be remarked that even this figurative usage has the erroneous 
connotation just mentioned. There are two other interrelated aspects of the general rela- 
tion ( |l"3D that are quite interesting and require further discussion. First, it is clear from 
equation (TT3|) that r + — t_ is nearly independent of Newton's gravitational constant G. 
Intuitively, this essentially comes about since in this approximation a "small" quantity 
is integrated over a "long" interval. Thus the result could be a "large" effect, since 
it is independent of the extremely weak gravitational coupling constant in this limit. 
It is useful to recall here a similar situation involving the total gravitational radiation 
energy emitted when a test particle of mass m that is at rest at infinity falls radially 
into a Schwarzschild black hole of mass M ^> m; the net result, ~ 10~ 2 (m / M)mc 2 , is 
also independent of G. Another analogous circumstance involves the net amplitude of 
relativistic nutation that is independent of c; that is, the spin vector of a test particle in 
an orbit of inclination a about a rotating mass undergoes - in addition to the normal 
precessional motions - a certain nutational motion of long (Fokker) period Tp ~ 2X0/3$ 
with an amplitude J sin a /Mr 2 uj that is independent of c JTB]. This is a consequence 
of the post-Schwarzschild approximation scheme that is discussed in section 4. 

The other significant aspect of the result ( |I3"D is that it is essentially independent of r for 
r ^> 2GM/c 2 . This gravitomagnetic clock effect is thus reminiscent of the topological 
Aharonov-Bohm effect. The connection between these effects can be further clarified as 
follows: The Aharonov-Bohm effect is simply related to the Sagnac effect via the Larmor 
theorem |17]]; in a similar way, the gravitational Larmor theorem provides a connection 
between the gravitomagnetic effect (|T^) and an analog of the Sagnac effect. This is 
discussed in the next section. 



3 Analogy with the Sagnac Effect 

A century ago, Larmor's theorem provided a local connection between magnetism and 
rotation. The theorem applies to circumstances involving slowly moving charged parti- 
cles and slowly varying fields whose strengths are considered only to first order. It turns 
out (cf. Appendix A) that the theorem can be extended to the gravitational case such 
that the gravitoelectric charge of a test mass m would be qp = —m and its gravito- 
magnetic charge would be qs = —2m. This is the content of the gravitational Larmor 
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theorem the negative signs of the gravitational charges account for the attraction 
of gravity and Qb/qe = 2 since gravity is a spin-2 field. 

It follows from Larmor's relation 

= (14) 

that in the gravitomagnetic case with qs = —2m, we have 

n L = -n P , (15) 

where ftp = B 9 /c is the local gravitomagnetic precession frequency of a gyroscope at 
rest. It follows from equation @ that in the equatorial plane (J-p = 0) one should 
obtain the same physical result at an orbit of radius p in the absence of gravity but in 
a frame rotating with uniform frequency ftp = GJ/c 2 p 3 . Imagine, therefore, two clocks 
moving in opposite directions with speed v on the circular orbit of radius p in the absence 
of gravity. According to an observer at rest in the frame rotating with frequency Ql, 
the periods of circular motion for the two clocks are 

t± = — 16 

by the nonrelativistic law for the addition of velocities. (The same result can be obtained 
with respect to static observers in the underlying inertial frame.) It follows that 



4(V)ii, 



which for v ^> pQi and v — *■ c reduces to the Sagnac effect for light - though the 
derivation presented here breaks down, of course, as it is only valid in the nonrelativistic 
approximation. That is, the Sagnac phase shift may be thought of - in the eikonal 
approximation - as the product of frequency of light and the total time difference for 
light to traverse the circular area in opposite directions. As is well known, the Sagnac 
effect is proportional to the enclosed area (up 2 ) ; in fact, the magnitude of the Sagnac 
phase shift is given in our case by 4uj*(7tp 2 )Ql/c 2 , where uo* is the frequency of the 
electromagnetic radiation. 

Let us now compute equation ( |17|) for the gravitomagnetic case v = puo, where u is 
the Keplerian frequency, u 2 = GM/p 3 , and 

| = ^ < 1 (18 > 

in all physically realistic situations. Hence we find that 

T+-T- ~ 47r4 = 47T-— , 19 

Uq Mc z 
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as expected from the gravitational Larmor theorem. 

The simplicity of the above argument is due to the fact that at a given p in the equatorial 
plane Ql is uniform; however, for an orbit at an inclination q^O this would no longer be 
the case. A more significant problem is that an orbit in the field of a rotating mass is not 
in general spatially closed. In fact, the equatorial circular geodesic orbits discussed thus 
far are exceptional in this respect. It is therefore necessary to formulate a generalization 
of relation fll3|) for an arbitrary orbit; this problem is treated in the next section. 



4 Standard Clocks in Space 

A spaceborne clock would in general follow a complicated orbit. Thus far, we have con- 
sidered only circular geodesic orbits in the equatorial plane of a rotating mass; however, 
it is possible to generalize the gravitomagnetic clock effect to the case of an arbitrary 
orbit. To clarify the situation, we first consider a "spherical" orbit that has a small 
inclination a< 1 with respect to the equatorial plane. Spherical orbits in Kerr space- 
time have been described by Wilkins ||19|| . Such a geodesic orbit is no longer spatially 
closed in general; therefore, it is necessary to define the relevant periods r± in terms of 
azimuthal closure. 

Most elementary astronomical systems, i.e. stars and planets, are nearly spherical 
bodies. The exterior gravitational field of a spherically symmetric mass distribution 
is uniquely described by the Schwarzschild spacetime; therefore, it is natural to ex- 
press the exterior field of such an astronomical system in terms of perturbations of the 
Schwarzschild field. The post-Schwarzschild approximation scheme is a useful method 



in relativistic celestial mechanics pig, |20| , gjfl . It is possible, in principle, to include 



the quadrupole and higher mass moments of the source in our treatment; however, the 
analysis of the gravitomagnetic clock effect for a general orbit would then become much 
more complicated and is beyond the scope of this paper. The unique first-order angu- 
lar momentum perturbation is simply given by the Kerr field linearized in the angular 
momentum parameter a. Thus we express the exterior field of a rotating mass by 

-ds 2 = -c 2 (l - 2<$>)dt 2 + (1 - 2<$>)- 1 dr 2 + r 2 (d6 2 + sin 2 9 dip 2 ) - 4cr^ sin 2 9 dt dip , 

(20) 

where $ = GJ/c 3 r 2 is a gravitomagnetic potential. For the exterior field of the Earth, 
$(r) < 7 x KT 10 and *(r) < 4 x 10~ 16 . 

Let us now consider a "spherical" orbit of small inclination « < 1 given by 

t = r (l-3^lW)r, (21) 
r = r-o, (22) 

9 = a sin r] + 3a — r (l — 2$ )co>rcos?7 , (23) 

2 c 

ip = (l - —T jLUT + ip , (24) 
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where $o = &( r o), k> is the proper frequency of the unperturbed orbit, u = ujqT , 

r = (l-3$ )-5 (25) 

and rj is - for the unperturbed orbit - the phase angle in the orbital plane measured 
from the line of nodes 

7] = UT + r] . (26) 

The constants (fo and 770 are given in terms of initial conditions at r = as in Figure 
1. The linear perturbation analysis is valid for u t <C (auo/c)^ 1 . It is possible to show 
that the average behavior of the orbit is described by the Lense-Thirring effect with a 
precession frequency VL LT ~ 2c^(r )/r . 




Figure 1: Plot of the direct "spherical" orbit of radius r . The azimuth of the line of 
nodes z/ is given by vq = (p a — arct an (cos at an 770). The polar angle at r = is given by 
cos#o — sinasin?7o. Here a is the inclination of the unperturbed (circular) orbit with 
respect to the equatorial plane of the rotating mass. 



It is clear from equations (^1|) - (0) that the orbit is no longer spatially closed once 
q^O. Let us define r + to be the period of motion from ip to (p + 2rr; hence, we obtain 
from equation (|24]) 



aw 

1 r )ujt + 

c 



2tt 



(27) 



It follows from this relation that (T = 27t/lv) 



t+ = T + 27rr c 



(28) 
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moreover, for a retrograde orbit, a — > —a so that 

t_ = T - 27rr - . (29) 

c 

Hence, we recover r + — r_ ~ Ana/c for $o <C 1. Thus the rotation of the source breaks 
the degeneracy in the proper period T of the orbit by an amount that is of the order of 
J/Mc 2 . This is illustrated in Figure 2. 




Figure 2: Schematic plot of clock orbits about a rotating mass. The azimuthal periods 
are illustrated for cos 770 > 0; in fact, the deviation in the polar angle after one orbit is 
A-naa cos 770/ cT . 



Let us next consider a general "spherical" orbit of arbitrary inclination a; in fact, the 
orbital equations have been derived in previous work pOj . We have 



t(r) = r ^1 — 3— r $o cosajr . 



r r 



r 



nf \ 1 ■ ■ \ 3au; _ , . .cosr? 

t/(r) = arccos(sm a sm 77) H 1 q( 1 — 2<Po) sm(2a) cut. 

2 c (7(77) 



acu 



<^(t) = arctan(cosatan77) H 

c 

where 77 is defined as in equation (|26l), (1(77) is given by 

a{rf) = (1 — sin 2 a sin 2 77) 



1 , , cos a 

2T 1 - 3r 1 - 2$ — T 



(30) 
(31) 
(32) 

cur + f , (33) 



(34) 
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and u is the longitude of the node as in Figure 1. For a€l, equations (^) - (|33|) reduce 
to the previous orbital equations - (|2"1|); moreover, the Lense-Thirring precession 



as well as other properties of the general orbit have been described in detail before 
For a = 7r/2, the main orbital motion is polar and the azimuthal motion indicates the 
Lense-Thirring precession of the nodes; therefore, the gravitomagnetic clock effect is 
absent in this case. Hence we assume a ^ tt/2 and consider r + from equation (|33|) such 
that ifio — > (fo + 2tt from r = to r = r + ; we find an implicit equation for r+ that can 
be solved perturbatively. Thus let 

r± = T ± 27r - A cos a ; (35) 

then, it is possible to show that 

A = T — 2Fq 1 tan 2 a cos 2 r] (36) 

for a sufficiently different from n/2 such that the perturbative treatment remains valid. 
It is important to recognize that A could change sign and become negative for sufficiently 
large inclination of the orbital plane with respect to the equatorial plane of the central 
source. For instance, for a and r/o both near 7r/4, A could go through zero and change 
sign. The azimuthal period of the prograde orbit would then become shorter than the 
retrograde orbit. 

The influence of the gravitomagnetic field of the source is in general reflected in any 
orbital timing mechanism. For instance, let us recall that t is the proper time of static 
observers at spatial infinity. Thus if the orbital motion about a rotating source is referred 
to static clocks at spatial infinity, we find 

t± = T ± 2rc - TqX' cos a , (37) 

c 



where T T = T and A' = A — 3 $ To using equations ( ffCD and (pop. This consideration 
could possibly be useful in the experimental determination of J/M as well. 

It follows from these considerations that in the approximation under consideration, where 
the rotation of the source is taken into account only to first order, the azimuthal orbital 
period r - as compared to the proper Keplerian period T - would be larger (smaller) 
due to orbital motion with A > in the same (opposite) sense as the rotation of the 
source. For the case of a satellite around the Earth, this difference is ±2na\ cosa/c, 
where 2na/c ~ 10~ 7 sec. A similar situation holds - as noted in the previous paragraph - 
even when the binary system is far away. It is interesting to contrast this circumstance 
with the variation in the period of the Hulse-Taylor binary pulsar PSR B1913+16. 
The binary period is monotonically decreasing by about 10~ 7 sec per orbit as a result 
of gravitational radiation damping, since the observations agree with the theoretically 
estimated decay due to the emission of gravitational radiation by the binary system. To 
measure an effect of this size for the motion of a clock in orbit around the Earth is the 
observational challenge posed by the gravitomagnetic clock effect (cf. section 6). 
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Though our approach has been macrophysical throughout, it is nevertheless interesting 
to note that for gravitational orbits around a neutron, say, r + — t_ is equal to the 
Compton period of the neutron. 

Finally, imagine standard clocks in orbit about a static astronomical source. The orbital 
period is defined by azimuthal closure as before. Suppose that for two clocks starting 
at r = in opposite directions around the source we have r + = r_. This degeneracy is 
removed once the source rotates [cf. equation ([35])]. It is then expected that r + — t_ = 
47raA cos a/c would be the dominant relativistic rotation-dependent term in general. 



5 Holonomy 

As already mentioned in section 0, the gravitomagnetic clock effect can be considered as a 
gravitational analog of the Aharonov-Bohm effect. It is well known that the Aharonov- 



Bohm effect is closely connected to the concept of holonomy (see, for example, [^2 
section 10.5.3.): The phase shift observed in the Aharonov-Bohm experiment can be 
obtained after the integration of [/(l)-parallel transport around a closed loop which 
surrounds a magnetic field. Mathematically, this phase shift is given by the holonomy 
attributed to a fibre bundle with base space S 1 (representing the closed loop), fibre U(l), 
and a [/(l)-valued connection (the vector potential that corresponds to the magnetic 
field). An analogous construction can also be conceived in order to obtain the gravit- 
omagnetic clock effect as a translational holonomy of a fibre bundle. In this case one 
defines a fibre bundle with base space S 1 (representing the parameter space of the az- 
imuthal coordinate ip), the real line R as fibre, and an i?-valued connection (describing 
the time shift of an orbiting clock due to gravitomagnetism). However, this construction 
is of limited use since the holonomy obtained is not directly related to gravitational 
holonomy. 

Gravitational holonomy is a characteristic of exterior parallel transport on spacetime. 
More precisely, it is attributed to an affine frame bundle which takes spacetime as base 
manifold and is equipped with an affine connection (cf. Appendix B). The gravitational 
holonomy itself splits into translational holonomy and rotational holonomy. In the limit 
of infinitesimally closed and contractible loops these holonomies turn, respectively, into 
the torsion and curvature of spacetime. In general relativity, one is restricted to Rieman- 
nian parallel transport such that torsion vanishes and curvature becomes Riemannian 
curvature. However, it should be noted that, in general, Riemannian parallel transport 
on a Riemannian manifold exhibits nonvanishing translational holonomy. This circum- 
stance does not imply the presence of nontrivial torsion, of course. Nevertheless, Petti 



[ f2~3|j employed spacelike circular curves in the Kerr geometry in a significant attempt 
to relate the angular momentum that occurs in the Kerr metric to torsion via transla- 



tional holonomy. Indeed, it was outlined in [23| how to assemble a manifold with torsion 



from torsion-less Kerr configurations in a manner similar to the assembly of a simplicial 
Riemannian manifold from piecewise flat spaces. This interesting undertaking deserves 
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further investigation since - as Petti [23] has observed - a careful limiting process should 
be taken into account. 

It follows from this that the gravitomagnetic clock effect does not give us any information 
on the local geometry of the Kerr metric that goes beyond Riemannian geometry This 
is not really surprising since an ideal standard clock, if considered as a structureless, i.e. 



pointlike, test mass, does not couple to torsion at all ||24|| ; in fact, it measures proper time 
along timelike worldlines and is thus intimately tied to the metric concept of Riemannian 
geometry. 



6 Discussion 

NASA's Gravity Probe -A involved the sub-orbital flight of a rocket carrying a hydrogen 
maser clock; the main result of this experiment was an accurate test of the gravitational 



redshift by Vessot et al. [[^, g6|. It is expected that some of the future space missions 
will carry ultra-stable clocks, and hence it is possible that the gravitomagnetic effect 
discussed in this paper could be measurable in the future. For an orbit of Keplerian 



period To, equations (35) and fl3"6|) imply that the relative gravitomagnetic variation in 
orbital period is 



T Mc 2 V P 3 
or about 4 x 1CT 11 for a near-Earth orbit. 

The principal aim of NASA's Gravity Probe -B will be the direct measurement of the 
gravitomagnetic field of the Earth by comparing the precession of gyroscopes with respect 
to telescopes on board a drag-free satellite in a polar orbit about the Earth jl3fl . After 



one orbital period, the gravitomagnetic precession angle of the gyroscope per 2% radians 
is 



n P 2GJ (GM 



u>o c 2 p 3 \ p 3 



-1/2 

(39) 

which is identical with equation (|3"8"p. This equality indicates that the theoretical possi- 
bility of measuring gravitomagnetism via clocks is essentially equivalent to the GP-B. 

The gravitomagnetic clock effect constitutes a surprisingly large effect. We mentioned 
already in section that in the case of an astronomical body such as the Earth, the time 
difference r + — t_ ~ AtcJ/Mc 2 after one orbit amounts to 

r + - r_ ~ 2 x l(T 7 sec . (40) 

Time shifts due to the gravitomagnetic field of the Earth are widely believed to be of 
much lower order. Indeed, suppose the time difference of a direct and retrograde moving 
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clock is taken at a fixed time, say, after one Kepler period To = 2k/ujq. Then it is 
straightforward to show that for a circular equatorial orbit, 



r + - t 




6 



(41) 



As an example, we set r = 7000 km and obtain 



t'+ -t'_~3x 10" 16 sec 



(42) 



This is nine orders of magnitude smaller than the time difference fllPf ). We recall that the 
result ( [40"D presupposes that the time difference of the two clocks is taken with respect 
to a fixed angle ip (i.e., after each clock has covered an azimuthal interval of 27r) and not 
with respect to a fixed time. 

In principle, it is a trivial task to measure a time difference of 2 x 10~ 7 sec with today's 
technology. However, an experimental verification of the gravitomagnetic clock effect 
does not only require the measurement of the time difference between two well-defined 
events up to an accuracy of 2 x 10~ 7 sec. It is, in fact, the proper time along the direct 
and retrograde orbits that is used as a clock. Therefore, it is not sufficient to send two 
highly accurate and stable clocks into space and let them orbit in opposite directions. 
The orbits themselves have to be highly accurate and stable as well. We recall that 
in order to obtain the time difference (fiGf) we have to subtract two periods, each of 
which represents the sum of a Kepler period and a much smaller relativistic correction. 
Under the assumption of identical orbits, the Kepler periods and their gravitoelectric 
corrections cancel upon subtraction while the gravitomagnetic contributions add up, 
yielding the actual clock effect under consideration here. Disturbances of the orbits will 
in general change the Kepler periods of the orbiting clocks. It follows that in this case 
the Kepler periods will not exactly cancel but may exhibit a significant time difference. 

In order to obtain some estimation of possible errors, we consider an actual experiment 
where atomic clocks are on board satellites in direct and retrograde orbits around the 
Earth. We may divide the error sources of such an experiment into two groups, namely 

(i) errors due to the tracking of the actual orbits, and 

(ii) deviations from idealized orbits due to 

— mass multipole moments of the Earth 

— radiation pressure 

— gravitational influence of the Moon, the Sun, and other planets 

— other systematic errors (e.g., atmospheric disturbances). 

The tracking of the actual orbits requires the measurement of distances and angles. 
The position of a satellite along an orbit can be determined to a few centimeters using 
the Global Positioning System (GPS); therefore, the temporal uncertainty that a near- 
Earth satellite has actually returned to the same azimuthal position in space is about 
5t ~ dr/v ~ 10 _6 sec. Here 5r « 1cm is the position uncertainty along track and v 
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is the orbital speed of the satellite. The gravitomagnetic clock effect, however, involves 
a definite temporal deviation of ~ 10~ 7 sec. It turns out that one should be able to 
measure the orbital radius up to an accuracy of the order of 10~ 2 cm and to determine 
angles up to an accuracy of 10 -10 rad in order to keep the errors due to the measurement 
smaller than the clock effect after one orbit. These requirements are about one order 
of magnitude higher than what can be achieved today. On the other hand, the clock 
effect is cumulative - just like the precession angle of a GP-B gyroscope - and hence 
many orbital periods can be used for a measurement of the gravitomagnetic effect; that 
is, the statistical tracking errors could be overcome if one were able to perform many 



single measurements [27 



The astrometric requirements for ensuring azimuthal closure should be emphasized. The 
operational definition of the azimuthal angle tp is ultimately based on the underlying 
astronomical coordinate system employed J^J. In fact, it may turn out to be necessary 



to monitor orbital motion so accurately as to be able to measure the Lense-Thirring 
effect and hence gravitomagnetism directly using only astrometric data on the orbits of 
spaceborne clocks. 

The systematic errors of the second group of error sources have a more serious influ- 
ence on the gravitomagnetic clock effect. In order to calculate the influence of such 
perturbative accelerations on the Kepler period one has to focus on sections of orbits 
rather than on complete closed orbits. This is because under favorable conditions (e.g., 
an almost constant perturbative acceleration) different temporal deviations can cancel 
each other when summed over a closed orbit. It is possible to estimate that perturba- 
tive accelerations should be kept below 10~ n g in order for the clock effect to become 
measurable. 

Perturbative accelerations due to multipole moments of the Earth are of the order of 
10 -3 g. The orbit of a satellite under the influence of the nonspherical and inhomoge- 
neous form of the Earth may be likened to a bumpy road. However, we do know the 
gravitational field of the Earth up to an accuracy of 10~ 9 g — 10 _10 g. NASA's already 
approved gravity mapping mission GRACE is expected to push this accuracy higher by 



about two orders of magnitude ||29|| . This would then make it possible, in principle, to 
correct for the influence of the multipole moments of the Earth on a gravitomagnetic 
clock experiment. 

The radiation pressure of the Sun causes perturbative accelerations of the order of 10~ 8 g. 
Using drag-free satellite techniques, this disturbance can be reduced by two orders of 
magnitude to 10~ 10 g. To keep this error source under control one thus has to be able 
to determine the solar radiation pressure accurately enough such that corrective calcu- 
lations can be performed. Alternatively, one must wait for drag-free satellites that can 
perform at least one order of magnitude better than current technology. 

The gravitational fields of the Moon and the Sun cause relative accelerations between 
the Earth and the orbiting clocks. The amplitudes of these accelerations are of the 
order of 10~ 7 g (Moon) and 10 _8 g (Sun). The influence of the other planets of the solar 
system plays only a minor role. For Jupiter, e.g., we obtain an influence of the order of 
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10~ 12 g. The positions of the Moon and the Sun are known with much higher accuracy 
than is needed to determine their gravitational field at the level of 10 -11 g; therefore, 
in principle, the influence of the Moon and the Sun can be properly taken into account 



It should be remembered that we have provided analytic expressions for t± only for equa- 
torial circular orbits and for spherical orbits with definite inclination. Our treatment 
must therefore be extended to eccentric orbits. This will probably require numerical 
integration of the corresponding geodesic equations. One should keep in mind that 
even initially circular orbits will acquire eccentricity over time due to the perturbative 
accelerations mentioned above. It is expected from these considerations and the previ- 
ous discussion that the measurement of this gravitomagnetic effect would be similarly 
complicated as in the GP-B. 

Finally it should be remarked that the possibility of sending electromagnetic signals 



around a rotating mass has been considered by a number of authors The time 

difference, according to the observer's clock, for a signal to follow a closed path (by 
means of "mirrors", etc.) in direct and retrograde directions around a rotating mass 
is proportional to the gravitomagnetic flux threading the loop. Similar effects arise in 
the synchronization gap around a rotating mass. One can show that all such effects 
are smaller than the clock effect considered in this paper by a factor of the order of 
$ = GM/c 2 r; $ is < 1CT 9 for the exterior field of the Earth and < 1CT 6 for the 
exterior field of the Sun [Q. The clock effect is larger than that for light by the inverse 
of a factor of the order v 2 /c 2 ~ $ ; intuitively, the clock accumulates a much larger 
gravitomagnetic effect since its Keplerian motion has a speed that is much smaller than 
the speed of light. 
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Appendix A: Gravitational Larmor Theorem 



The local Larmor equivalence between magnetism and rotation has an exact analog in 
the theory of gravitation. To elucidate this connection, let us consider the motion of free 
test particles in the gravitational field given by the metric form (1). The geodesic motion 
follows from 5 J —mc ds = ; however, for purposes of comparison with the Newtonian 
equations of motion j32|, we consider 5 J Ldt = with L = —mcds/dt . That is 

1/2 



-m c 



2$ 



2$)- + -vA 9 



(Al) 
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It follows from this Lagrangian that the canonical momentum is given by 



P = mV(l + 2$) 



~dt 



2m 



c 



TA. 



(A2) 



where T = cdt/ds and r w 1 for nearly Keplerian motion under consideration here. The 
analog of equation (Al) in electrodynamics is 



where <fi is the electric potential. A perturbative treatment of equation (Al) can bring it 
to the form of equation (A3) with q — > —m except for an extra factor of 2 in the term in- 
volving the gravitomagnetic vector potential; indeed, the origin of this extra factor is that 
off-diagonal terms appear in the spacetime interval with a factor of 2. The simplest way 
to deal with this situation is to assign a gravitoelectric charge of — m and a gravitomag- 
netic charge of —2m to the test particle. In this way, equation (A2) becomes consistent 
with the standard relationship in electrodynamics between the canonical momentum P 
and the kinetic momentum p of a particle of charge q in an inertial frame in Minkowski 
spacetime, P = p + gA/c. It is convenient to use units such that G = 1 in what follows. 
In these units, let us therefore assign gravitomagnetic charges of qp = —2m to the test 
particle and Qp = 2 M to the central source. The source is rotating; therefore, it has 
a gravitomagnetic dipole moment fi g = Qp3/2Mc = 3/c and hence a gravitomagnetic 
vector potential A g = fi g x pj p 3 . This is consistent with equation (2); hence, the 
analogy established thus far between magnetism and gravitomagnetism is exact in the 
approximation under consideration here. That is, the gravitomagnetic dipole moment 
for a test gyro of spin S is then p, g = — S/c and the rate of precession of this dipole 
moment in the gravitomagnetic field, dS/dt = p, g x B g , is precisely given by equation 
(§). Moreover, the comparison of equation (Al) with its electrodynamic analog (A3) 
leads to the interpretation that the factor of 2 in the gravitomagnetic charge has to do 
with the fact that ds 2 is a quadratic form or, equivalently, that gravitation is a spin-2 



The correspondence between magnetism and gravitomagnetism can be simply extended 
to the Larmor theorem. Thus the gravitomagnetic Larmor frequency is given by dp — 
q~Q g /2mc = —~Q g /c = —dp . The gravitational Larmor theorem obviously holds for the 
gravitomagnetic precession of a gyroscope, since it means that the same rate of precession 
would be obtained locally for a free ideal gyroscope in the absence of gravitation but 
observed from a frame rotating with frequency flp = — ftp. Further applications of 
these ideas are contained in fT8| . 

In view of the intrinsic significance of the gravitational Larmor theorem and its basic 
relationship with Einstein's principle of equivalence, it is important to point out some 
of its underlying features. The theorem is formulated with respect to preferred standard 
observers in the background Cartesian coordinate system, i.e. the rest frame of the center 
of mass of the rotating source. Furthermore, it has been assumed in our discussion thus 
far that gravitomagnetism arises from mass currents, as in the case of a rotating planet 




(A3) 



field. 



16 



or star, so that the gyromagnetic ratio is then unity. In this case, qE ^ 1b and hence 
Lorentz invariance is broken. It is necessary to point out, however, that our treatment 
can be extended to the gravitomagnetic field of a rotating black hole as well. 

Consider first a test gyro in the field of an arbitrary rotating configuration. It is reason- 
able to suppose that the gravitomagnetic field is independent of the specific nature of 
the source in the linear approximation under consideration here. Hence, let \x g = kJ/c 
for the source and p, g = —kS/c for the test gyro, where k is a universal constant re- 
gardless of the nature of the rotating systems. Then, dS/dt = fl g x B 3 implies that 
k 2 = 1 and we choose k = 1 by convention and without any loss in generality. Next, let 
/i = r yQJ/2Mc, where 7 is the gyromagnetic ratio and k = ^Qb/ZM =1. In view of 
our previous arguments, cf. equations (Al)-(A3), Qb = 2M and 7 = 1; however, in the 
absence of such arguments the simplest possibility is Qe = Qb = M and hence 7 = 2. 
In fact, the charged Kerr system has /i = Qa and 7 = 2. In this way, the gravitational 
Larmor theorem can be extended to the exterior field of a Kerr system. Our results sug- 
gest, but do not prove, that 7 = 2 for the Kerr- Newman system since gravity is a spin-2 
field. That is, the factor of 2 in the gravitomagnetic charge of a rotating mass that is 
connected with the spin-2 character of the linearized gravitational field simply goes over 
to the gyromagnetic ratio for the case of a charged Kerr configuration. This apparent 
independence of 7 from the detailed electromagnetic nature of the source for this case 
is supported by the circumstance that 7 = 2 extends to the generalized Kerr-Newman 
spacetime |j33| |, which contains an infinite set of multipole moments, indicating the in- 



sensitivity of the gyromagnetic ratio to the detailed structure of the source. Further 
evidence for this viewpoint comes from the fact that a proper classical interior electro- 
magnetic source for the exterior Kerr-Newman field should have 7 = 1 in the absence of 
gravity. In fact, the gyromagnetic ratio of the Kerr-Newman system loses its meaning 
once one isolates the electromagnetic source of the Kerr-Newman geometry by simply 
turning the gravitational interaction off; that is, one is then left with only a charge Q 
and a static magnetic dipole moment Qa. 



Appendix B: Gravitational Holonomy 

In general relativity the local geometry of spacetime is determined by a metric g. Parallel 
transport in spacetime is then described by means of a Riemannian connection one-form 

r oP which is derived from the metric g. The prescription of parallel transport is the 
following: Consider a linear frame e Q at a point x of spacetime. Parallel transport of 
e a (x) to an infinitesimally neighboring point x' = x + dx is defined by the identification 
of the frame 

(e a + de a )(x') = e' a + de a (x') := e' a + Tj(x')e'g (Bl) 
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with the parallel transported frame e a (a;).[] It is well known that this Riemannian parallel 
transport is both metric preserving and torsion-free. The Riemannian curvature two- 
form 

Rj := dT a p - IV A r/ (B2) 

is a local measure of the nonintegrability of Riemannian parallel transport. That is, if 
we parallel transport a frame around an infinitesimal closed loop in spacetime the frame 

{} . 

will undergo a "pseudo-orthogonal" linear transformation which is determined by R a . 
We can generalize this concept and consider parallel transport of a frame around a finite 
loop. This will again result in a "pseudo-orthogonal" linear transformation of the frame 
and constitutes the holonomy. More exactly, it is the holonomy of a linear frame bundle 
with spacetime as base manifold and a Riemannian connection as connection. 

More general non- Riemannian spacetime geometries are determined not only by a metric 
g but also by an independent affine connection which usually can be taken as a Cartan 
connection {^"^Yj 5 ) [34]. The Cartan connection determines affine parallel transport 



which generalizes Riemannian parallel transport. While Riemannian parallel transport 
determines the identification of neighboring linear frames by means of (pseudo-) ortho- 
gonal transformations, the affine parallel transport determines the identification of neigh- 



boring affine frames by means of affine transformations [35|. Note that an affine frame 
at a point x of a spacetime M is a pair (e a ,p)(x) which consists of a linear frame e a (x) 
and a point p(x). The affine frame is defined as an element of the affine tangent space 
A X M. We recover from the affine frame (e a ,p)(x) the linear frame e a (x) if and only if 
p E A X M is identified with x. 

The precise prescription of affine parallel transport is the following: Consider a linear 
frame e a at a point x of spacetime and define as := i? a — 8fdx % . Affine parallel 
transport of (e a ,p)(x) to an infinitesimally neighboring point x' = x + dx is defined by 
the identification of the affine frame 

(e Q + de a )(x') = e' a + de a (x') := e' a + T P{x')e' p E A X >M (B3) 

(p + dp){x') = p + dp(x') := p' + r (T)c V) e' a E A X ,M (B4) 

with the parallel transported frame (e a ,p)(x). Integration of (B3) and (B4) around an 
infinitesimal closed loop yields the torsion and curvature two-forms T a and R a . The 
explicit defining formulas are 

torsion T a := Dti a = d$ a + Tp* A , (B5) 

curvature Rj := dT a p - T a 7 A T 7 /3 . (B6) 

Torsion and curvature determine the resulting affine transformation of an affine frame 
after parallel transport around an infinitesimal closed loop. The torsion measures the 



1 This definiton does depend on the choice of e' a — e a (x'), of course. Another choice of e' a would 
correspond to another gauge of r J 3 . 
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translational part of this transformation while the curvature measures the homogeneous 
part. Generalizing affine parallel transport to parallel transport around finite loops 
yields the translational and rotational holonomy which we referred to as gravitational 
holonomy. It is the holonomy related to an affine frame bundle with spacetime as base 
manifold and a Cartan connection as connection. Further discussion of gravitational 
holonomy can be found, for example, in references p6l Wj\ . 
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